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Abstract 

Discrete flavour symmetries have been proven successful in explaining the lcptonic flavour 
structure. To account for the observed mixing pattern, the flavour symmetry has to be 
broken to different subgroups in the charged and neutral lepton sector. However, cross- 
couplings via non-trivial contractions in the scalar potential force the group to break to 
the same subgroup. We present a solution to this problem by extending the flavour group 
in such a way that it preserves the flavour structure, but leads to an 'accidental' symmetry 
in the flavon potential. 

We have searched for symmetry groups up to order fOOO, which forbid all danger- 
ous cross-couplings and extend one of the interesting groups A4, T 7 , S4, T' or A(27). 
We have found a number of candidate groups and present a model based on one of the 
smallest extensions of A4, namely Qs X A4. We show that the most general nonsuper- 
symmetric potential allows for the correct vacuum alignment. We investigate the effects 
of higher dimensional operators on the vacuum configuration and mixing angles, and 
give a see-saw-like UV completion. Finally, we discuss the supersymmctrization of the 
model. Additionally, we release the Mathematica package Discrete providing various 
useful tools for model building such as easily calculating invariants of discrete groups and 
flavon potentials. 
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1 Introduction 



Over the last decade, neutrino oscillation experiments have measured the mixing angles of 
the leptonic mixing matrix Upmns to quite some accuracy [1-3]. It turns out that two 
of the mixing angles, namely the solar and atmospheric angles #12 and 023, are large while 
the third one, the reactor angle #13, is small. Recently, there has been a hint of a non- 
vanishing third mixing angle #13 by the T2K experiment [4] close to the upper bound of the 
CHOOZ experiment [5]. Remarkably, the current best fit values (taking into account the 
recent measurements) for the case of normal neutrino mass hierarchy [1] 

sin 2 0i2 = 0.312+°;°^, sin 2 023 = 0.52+°;°? and sin 2 6 13 = 0.013±g;gg| 

are rather close to tri-bimaximal mixing 

sin 2 0i2 = -, sin 2 023 = - and sin 2 0i3 = , 

3 2 

which was first proposed by Harrison, Perkins and Scott [6] . Before the hint of a non- vanishing 
013 mixing angle by T2K and MINOS [ ], the description by tri-bimaximal mixing was even 
better. 

This interesting observation has led many authors to search for a theoretical explanation 
of this curious mixing pattern. While in the quark sector the mixing angles are small and 
can be explained by a Froggatt-Nielsen type U(l) symmetry [8], which also accounts for the 
quark mass hierarchies, in the neutrino sector so far the most fruitful approach has been to 
introduce a non-abelian discrete symmetry, like A4 [9; 10], T7 [11], £4 [12; 13], T' [14] and 
A(27) [ ] among others. The left-handed lepton doublets are commonly assigned to a non- 
trivial representation of the flavour group, e.g. a triplet. Subsequently, this symmetry is 
spontaneously broken to different subgroups in the charged lepton and neutrino sector [16; 17]. 
This requires at least two scalar fields to obtain vacuum expectation values (VEVs), which are 
pointing in two different directions in the space of the flavour symmetry. This is commonly 
denoted by VEV alignment, which we are going to address in this article. 

The most well-studied non-abelian discrete group is the group A4, which is the symmetry 
group of a regular tetrahedron. It is the smallest group with an irreducible three dimensional 
representation 3. If one assigns the lepton doublets as L ~ 3 and charged leptons to the 
three singlet representations e c , /U c , r c as I2 anc ^ l3i the tri-bimaximal mixing structure 
is generated due to a mismatch of the vacuum expectation values of the flavons x ~ 3 an d 
4> ~ 3 that couple to charged leptons and the neutrinos, respectively. If A4 is broken down 
to the subgroup Z3 in the charged lepton sector, e.g. by (x) oc (1, 1, 1) T , and it is broken to 
the subgroup Z2 in the neutrino sector, e.g. by (</>) oc (1,0, 0) T , the unitary transformation 
that connects the most general mass matrices invariant under these symmetries is given by 



Upmns = U H psU2z{9) with U H PS = —j= —j= 

v -L 4* -L 1 

\ V6 V3 V2 / 



(1) 



2 



where Uhps denotes the tri-bimaximal mixing matrix and U2s(9) denotes a 2 — 3 rotation by 
the angle 9, which is generated by operators of the form LL{yq — h.c.) with r/ ~ !2,3- I n 
case, these operators do not contribute to the neutrino mass matrix, i.e. Im (y (rj)) = 0, the 
leptonic mixing matrix is given by the tri-bimaximal mixing matrix Uhps- 0- n light of the 
recent hint on a non- vanishing #13, a non- vanishing operator of this type has been discussed 
in [18].) 

The prediction of tri-bimaximal mixing in A4 models thus requires a special vacuum 
alignment 1 , which should have a dynamical origin within the model. In the most straightfor- 
ward dynamical model, namely the usual scalar potential, the cross coupling terms connect- 
ing <fi and x Via non-trivial A4 contractions, e.g. (x 2 )l (^ 2 )l Q > forbid the desired vacuum 

— o 

alignment [16; 17], as will be reviewed in the next section. This vacuum alignment problem 
is not limited to models, but a general problem of most of the symmetry groups, that 
have been studied. In the literature, several mechanisms have been proposed to address this 
problem: 

(i) in models with extra dimensions, it is possible to localize the two flavons differently 
in the extra dimensions with no (or negligible) overlap of the wave functions, thereby 
forbidding (or suppressing) all cross-couplings (see e.g. [10] for a localization on different 
branes) 2 , 

(ii) in a supersymmetric framework it is possible to introduce R charges and driving fields 
with R charge 2 such that the tadpole of the driving fields forces the desired vacuum 
alignment (see e.g. [16]). 

While these approaches are interesting and worth studying, we will develop further an idea 

by Babu and Gabriel [21], who have suggested a group-theoretical mechanism to forbid the 

dangerous cross-couplings and thus a way to realise the VEV alignment without using R- 

symmetries in supersymmetry or brane constructions. 

They proposed an extension of the flavour group A4 in such a way that the Standard 

Model leptons only transform under the A4 subgroup of the full flavour group. In the scalar 

sector, the flavon x °f the charged lepton sector also transforms only under the A\ subgroup, 

while the flavon <f> of the neutrino sector transforms under the full flavour group G. For a 

suitably chosen group G, it is then possible that the additional group transformations forbid 

the contractions {(f>(j>)\ and (<j)4>)^ , which lead to the dangerous couplings in the scalar 

2,3 1 

potential and make the correct vacuum alignment impossible. In other words, the additional 
discrete symmetry leads to an accidental symmetry at the renormalizable level in the flavon 
potential, G x A4, allowing for a different breaking of the two A4 subgroups of the accidental 
symmetry. The coupling to leptons only respects the diagonal A^ subgroup, which is thus 
broken to different subgroups in the charged lepton and neutrino sectors, as desired. 

1 See [19] for a model that can accommodate the large value for #13 suggested by T2K and still needs a 
special vacuum alignment. 

2 Models with extra dimensions also allow an explicit breaking of the flavour symmetry via boundary con- 
ditions [20]. 
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Note that this construction requires that the additional group generators cannot all com- 
mute with the generators of A±, i.e. the flavour group G cannot be a direct product of A4 
with some other group, but it has to be a slightly more general object, a semidirect product. 
In fact, we will generalize this construction one step further and look into more general group 
extensions. This will be explained in more detail in sec. 3. 

In their work, Babu and Gabriel used a special type of semidirect product, a so-called 
wreath product of A4 with 53, i.e. the product of four factors of £3 which are evenly permuted 
by the group A4. It is thus a very complicated flavour group of order 12 • 6 4 = 15552 and 
requires the use of very large representations up to dimension 48. This model further suffers 
from a fine-tuning problem, as the diagonal and off-diagonal elements of the neutrino mass 
matrix are generated by operators with very different mass dimension, 5 and 10, even though 
both entries should be of comparable size. 

In this article, we are addressing these issues and are presenting the result of a search 
for simpler and more attractive semidirect product groups G = N x H as well as general 
group extensions G satisfying G/N = H with H being A4, T7, S4, T" or A(27) 3 , which 
lead to an accidental symmetry G x H in the flavon potential at the renormalizable level. 
We included all discrete groups up to order 1000 in our search and found several candidate 
groups. The smallest candidate groups are of order 96, in particular the semidirect product 
group of the quaternion group with A4, Qs x A4, which we discuss in more detail in sec. 4. 
This group does not have representations of size larger than four and, in the model we present, 
on-and-off-diagonal entries in the neutrino mass matrix are generated at the same order. 

This work has been accompanied by the development of the Mathematica package Discrete 
facilitating the calculation of the different covariants of a discrete flavour group. Dirac and 
Majorana mass matrices as well as the flavon potential can be calculated automatically up 
to an arbitrary order. It can access the large group catalogues implemented in GAP [ ]. Its 
SmallGroups [ 5] catalogue, for example, contains all discrete groups up to order 2000 with 
the exception of order 1024. 

The outline of the paper is as follows. In sec. 2, we discuss the VEV alignment problem 
in the context of A4. Our search for semidirect product groups and the results are described 
in sec. 3. Readers, who are not interested in the technical details of the construction, may 
skip this section. In sec. 4, we discuss the smallest candidate group Qs x A4 and construct 
a model using it. It is the minimal model based on A4 allowing for the correct vacuum 
alignment. Higher order corrections are discussed in sec. 5. An ultraviolet (UV) completion 
is presented in sec. 6 and a supersymmetric version is given in sec. 7. The Mathematica 
package Discrete is introduced in sec. 8. Finally, we conclude in sec. 9. Group-theoretical 
details are summarised in the appendix. 



3 This mechanism is of course not limited to these five groups, but is also relevant for other flavour groups, 
e.g. larger groups constructed to be in agreement with the recent T2K measurement such as [19]. 
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(b) Generators 
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Table 1: Character table of A4 as well as matrix representation of generators in 3. Here 
oj = e i2n / 3 . 

2 VEV Alignment in A 4 Revisited 

In this section, we want to remind the reader about the difficulties one encounters when 
minimising flavon potentials [16; 17]. Here we focus on the problem one faces in the most 
straightforward case - namely the case of a nonsupersymmetric scalar potential. The case of 
softly-broken supersymmetry is included in this analysis as SUSY only further restricts the 
dimensionless couplings of the potential while care has to be taken not to have flat directions 
in the cubic superpotential. We will come back to the SUSY case in sec. 7. 

For simplicity, we consider A4, the symmetry group of the tetrahedron, which is the 
smallest discrete group with a three dimensional irreducible representation. It is presented by 
(S, T\S 2 = T 3 = (ST) 3 = l). As we have discussed in the introduction, tri-bimaximal mixing 
is generated by breaking this group to its subgroups generated by S and T in the neutrino 
and charged lepton sectors, respectively. The character table and the representation matrices 
for the three dimensional representation are given in Table 1. 

Let us look at the potential 4 



Vy 



m l {xx)l x + Ai (xx)l 1 (xx)l 1 + A 2 (xx)l 2 (xx)l 3 



(2) 



of a real scalar triplet % °f ^4 that couples to charged leptons via the operators £Hf c x/A- and 
should therefore acquire a VEV (x) = (v , v , v ) conserving the Z% subgroup generated by 
T. The symmetry breaking of A4 to {S) = Z2 (which might or might not be due to the VEV 



of another triplet <ft in the neutrino sector with ((f)) 
soft terms in the potential: 



(w, 0, 0) ), will lead to the following 



soft,Z 2 



2 2, 2 2, 2 

m A Xi + m B x 2 + m c X2X3 



(3) 



The minimisation conditions of the full potential V = V x + V so ft,z 2 evaluated at the desired 



4 Here, we have assumed a discrete symmetry (\ -> — Xi f c ~~f c > with / = e,fi,r) that separates the 
charged lepton from the neutrino sector, as it is common practice. The operator (xx)s^ ' (xx)s^j which one 
would naively expect, can be expressed as a linear combination of the other operators. 
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minimum (%) = (v', v' , v') result in 
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ml + V3m 2 A ^j v' + 4X lV ' 3 



V 



dX3 
d 



2m| v' 



9xi 9x3 



V 



Xi=v 



Xi=V 



(2 m\ — ttiq) v 



(4a) 
(4b) 
(4c) 



The vacuum alignment thus requires m 2 B = and m 2 c = 2m\ and therefore two completely 
different A4 contractions need to have the same coupling in the scalar potential, an option we 
exclude as fine-tuning. Even if one sets the terms m\ B c to zero, they will still be generated 
on loop-level and disturb the VEV alignment. The breaking of A4 to two different subgroups 
thus requires a systematic mechanism to forbid m 2 A B c . 

On the contrary, soft breaking terms which preserve the same Z3 subgroup 



V f 



soft,Z 3 



(X2X3 + X1X2 +X3Xl) 



(5) 



are not in conflict with the VEV alignment, because they do not change the structure of the 
minimisation conditions at the minimum (%) = (v, v, v) 







9V X + Vsoft,z 3 
9xi 



\v (3m 2 s + V3m 2 + 6Aii/ 2 ) 



(6) 



Hence, the flavon potential enforces the VEVs to align. 

We thus conclude that the desired vacuum alignment requires a mechanism to forbid 
all couplings between the flavon sectors that break A4 to Z2 and Z3, respectively, except 
for the quartic coupling where both couple in pairs to singlets. This can be rephrased in 
the requirement to have an 'accidental' symmetry A\ x A^ in the flavon potential, where 
the first group factor, A%, corresponds to the flavons coupling to neutrinos and the second 
one, A4, to flavons coupling to charged leptons. Note, that the Kronecker product 3x3 = 
II + 12 + l3 + %s + 3 A allows couplings of the form (xx)l 2 (#)l 3 , (xx)l 3 (#)l 2 an d 
(xx)s( c t >c l ) )s m the minimal A4 model discussed above and the desired vacuum alignment is 
thus not possible. These kind of couplings can not be forbidden by assigning % or <f> to a 
unitary representation of an additional internal symmetry group commuting with the flavour 
group, because x^X an d <$ ( t ) wm always be invariant. In particular, it is not possible to solve it 
by introducing an additional commuting group factor, which is a discrete group or a compact 
Lie group. 

For this reason, the VEV alignment problem has been mainly studied within the context 
of SUSY as well as brane constructions within extra-dimensional models circumventing this 
problem. In the following section, we show how the required vacuum alignment can be 
achieved with an internal symmetry group by extending the flavour group in a non-trivial 
way. 
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3 Group Extensions and Vacuum Alignment 



In the following, we explain the type of groups we are searching for and why we are searching 
for these groups. We always use the group as an example, but the arguments hold for any 
group. In a first step, we directly extend the group by adding new generators, which do not 
commute with the generators of the flavour group. In the second subsection, we generalize 
our approach and look for general group extensions. The group theoretical notions we use 
are defined in the footnotes of this section. 

3.1 Semidirect Product Groups 

To reproduce the success of A4 models, we search for an extended flavour group 

G = (S,T,X U . . .,X n \S 2 = T 3 = (ST) 3 = r*(Xi, ...,*„) = rf x (S, T, X u . . . , X n ) = l) 

(7) 

that contains H = (S, T\S 2 = T 3 = (ST) 3 = l) ~ A4 as a subgroup. Xi denote the additional 
generators of the extended group and , r^ lx with a = 1, . . . , s a and (3 = 1, . . . , sp additional 
relations. Note that there are no additional relations involving only S and T. As we discussed 
in the last section, not all of the additional generators can commute with H. Therefore, there 
have to be relations r™ lx . These generators will be needed to forbid the dangerous couplings 
discussed in the last section. Any discrete group that contains H as subgroup can be written 
in this way. 

We further demand that there should be representations p\ : G — > GL(V) with 

Pi (X 3 ) = t Vj = l,...,n (8) 
and pi(S) and p\(T) corresponding to the usual A± representations i = 1±, I2, I3 and 3, e.g. 



( 1 

0-1 I , ps(T) 



0-1 



/ 1 

1 I . (9) 
^10 



If the SM fermions are assigned to these representations, the A4 predictions for the mixing 
angles remain unchanged. The existence of the representation p = p% gives a first constraint 
on the flavour group G: The image of the representation p is isomorphic to H, i.e. im(p) = H, 
and its kernel 5 = ker p = {Xi, ...,X n ) is a normal subgroup 6 of G with the quotient group 7 
G/N = im(p) = H (by the first isomorphism theorem), p thus essentially defines a surjective 
homomorphism 8 from G onto H, which is the identity on H and whose kernel is N. Groups of 



5 The kernel of a representation p is defined by kerp = {g £ G\p(g) = 1}. 

6 A normal subgroup TV of a group G, denoted by N<\ G, is a subgroup, which is invariant under conjugation 
by an arbitrary group element of G, i.e. gNg^ 1 — TV. 

7 The quotient group G/N is defined by the set of the left cosets gN with g G G. 

8 A (group) homomorphism p : G — > H is a mapping preserving the group structure, i.e. p(gi<?2) = 
p(ffi)p(S2) V<?i,2 G G. A surjective homomorphism p : G — > H has the additional property im(p) — H. 
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this type are known as semidirect product groups 9 G = N x H, which is a generalisation of 
a direct product N x H. One example of a semidirect product group is A4 = [Zi x Z2) x Z3 
itself. As N and H can not commute, G can not be a direct product, G ^ N x H . 

Once we have found such a group we can assign the lepton doublets, charged leptons and 
the flavon % that couples to the charged lepton sector in the usual way to representations 
3 and lj, i = 1,2,3, while assigning the navon <p of the neutrino sector to an irreducible 
representation of G, which is faithful 10 on N, and contains 3 in the Kronecker product 4> n at 
some order n. If this representation <f> was not faithful on N, it would be possible to restrict 
to the smaller group G/ker0|jv (by the third isomorphism theorem), which leads to the same 
flavour structure, and study its predictions. The problematic cross-couplings (xx)i2(00)lg> 
(xx)l (<W>)l an d (xx)3(00)3 can now be forbidden, provided that the Kronecker product 
<p x <fi does not contain the representations 3 as well as X2,3- Thus, the navon potential of (f> 
and x exhibits an 'accidental' symmetry G x H at the renormalizable level. This accidental 
symmetry is broken to G at an higher order in the navon potential. 

We thus systematically search for flavour groups G containing a subgroup H and a normal 
subgroup iV satisfying G/N = H{= A4), which lead to an 'accidental' symmetry G x H in 
the renormalizable part of the navon potential. Using the computer algebra system GAP [22] 
and its SmallGroups catalogue [23], we have performed a scan over all discrete groups G 
up to order 1000. As the vacuum alignment problem is not specific to the group A4, we 
have searched for semidirect product groups N xi H with the desired properties for the groups 
H = A4, T7, S4, T' and A(27), which are known to be interesting for flavour model building 11 . 
We applied the following conditions: 

1. G=NxH^NxH with H being one of the groups A 4 , T 7 , 5 4 , T or A(27); 

2. there is an irreducible representation (ft, which is faithful on N; 

3. (j) n contains 3 for some n; 

4. there is an 'accidental' symmetry G x H in the renormalizable part of the navon po- 
tential, i.e. there are only couplings via the trivial singlet between x an d <ft a t the 
renormalizable level, e.g. only {x 2 )\^ ( ( / )2 )l-^ exists for real representations x> 4>\ 

It turns out that there are only candidates for A4, T' or S4 up to order 1000, which are 
presented in Tab. 2. Although, there are semidirect product groups which fulfil the first three 
criteria for H = Tj, or H = A(27), none of them leads to the desired accidental symmetry in 

9 A group G is a semidirect product of a subgroup H and normal subgroup iV if there exists a homomorphism 
G — > H which is the identity on H and whose kernel is N. The direct product N x H can be considered as a 
semidirect product, where H is a normal subgroup of G as well. The two factors N and H of a direct product 
commute. 

10 A representation cj> is faithful, if the homomorphism <f> : G — > GL(V) is injective. It is faithful on a 
subgroup N, if 4>\n is faithful. 

n All of these groups have size less than 30 and contain a three-dimensional representation. There are two 
more groups with this property: A4 x Z2 and Z g x Z3 [2 I]. See also [25] for a recent overview of finite groups, 
which are useful for flavour model building. 
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Subgroup H 


Order of G 


GAP 


Structure Description 


Z{G) 




96 


204 


Q 8 x Ai 


z 2 




288 


860 


T x A 4 


z 2 




384 


617, 20123 


((Z 2 x Q 8 ) x Zo) x Aa 


z 2 




576 


8273 


{Z 2 .S 4 ) X A4 


z 2 




768 


1083945 


(Z 4 .Z|) x A 4 


z 4 




1085279 


((Z 2 x Qi 6 ) x Z 2 ) x A 4 


z 2 




192 


1494 


Q 8 x S4 


z 2 




384 


18133, 20092 


(Z 2 x Q 8 ) x S 4 


z 2 






20096 


((Z 4 x Z 2 ) x Z 2 ) x 54 


z± 




576 


8282 


T' x S 4 


z 2 


5 4 


8480 


(Z 3 x Q 8 ) x S 4 


z 6 




768 


1086052, 1086053 


((Z 2 x Q 8 ) x Z 2 ) x S 4 


z 2 




960 


11114 


(Z 5 x Q a ) x 5 4 


Zio 




192 


1022 


Qs x T' 




V 


648 


533 


A(27) x T' 


Z3 


768 


1083573, 1085187 


((Z 2 x Q 8 ) x Z 2 ) x T 
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Table 2: Candidate groups G up to order 1000 that may be written as non-trivial semidirect 
products G = N x H for the groups H = T7, S4, T' , A (27) and that lead to an enhanced 
symmetry in the scalar potential making the correct vacuum alignment possible. No such 
groups were found for H = T7, A (27). Details of the groups may be accessed using the 
computer algebra system GAP by using the command SmallGroup( Order, GAP). Qs denotes 
the quaternion group, which is defined in sec. 4 and the generalized quaternion group of order 
16, Qi6, is defined by Q±q = (x,y\x 8 = 1, x 2 = y 4 , y~ l xy = y^ 1 ). The expression of the form 
N.H is the GAP notation of a central extension, i.e. N is a normal subgroup of G, which is 
contained in the centre of G, and H is the quotient group G/N = H. Note that there can be 
more than one semidirect product of N by H. 

the scalar potential. This might be related to the fact that these groups have complex three- 
dimensional representations, and there are more couplings that would have to be forbidden 
by the additional symmetries than in the case of H = A4, T" and S4, which have real three 
dimensional representations. Additionally, there are simply less groups up to order 1000, 
which can be considered as an extension of T7 or A(27) compared to the other groups. 

Looking at the list of candidate groups, we further note that the normal subgroup N 
is non-abelian for all our candidate groups. In addition, the defining homomorphism 12 of 
each semidirect product is injective for H = A 4 , S4 3 and in case of H = T 1 , each group 
iV x X" allows for a defining homomorphism with image A4 or T". The quaternion group Qs, 
which frequently appears in Tab. 2, is the smallest non-abelian group allowing for a defining 

12 Equivalently to the previous definition, a semidirect product N x H can be defined via a homomorphism 
if : H — > Aut(iV), where Aut(iV) denotes the group of all automorphisms of N, i.e. the isomorphisms iV — > TV. 
The defining homomorphism is sometimes indicated as index of x, i.e. N x v H. 

13 The same applies for the wreath product S3 x A4 introduced by Babu and Gabriel [21]. 
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homomorphism with these properties. Furthermore, all candidate groups have a non-trivial 
centre 1 Z(G) <3 N. Hence, the representations can be classified according to their way of 
representing the elements in the centre, i.e. whether (a subgroup of) the centre is represented 
trivially (mapped to the identity) or not. In particular, the representations x °f G, which 
are directly related to irreducible representations x °f H with x\h = X H map the centre to 
the identity. They are single valued (in analogy to the representations of SU(2) with integer 
spin). However, groups that fulfil these conditions do not necessarily have to have a non-trivial 
centre. For example the wreath product S3 x A4, introduced by Babu and Gabriel [21], has 
a trivial centre. 

Before studying the vacuum alignment for the smallest candidate group, let us look more 
closely at how the breaking to different subgroups leads to the flavour structure. As has been 
mentioned in the introduction, it has been argued in [16; 17] that the neutrino mixing matrix 
can be obtained by breaking the flavour group to different subgroups in the charged and 
neutral fermion sector, respectively 15 . It is usually broken by flavoured scalar fields acquiring 
a VEV, which breaks the flavour group to the corresponding little group 16 . However, the little 
group of the VEV of a scalar field is not necessarily the little group relevant for the flavour 
structure, as the mass term might be generated by a Kronecker product of several scalar 
fields, e.g. the neutrino mass matrix might be given by £HIH ((ft) 2 and, therefore, the little 
group of ((ft) 2 is the relevant one 17 . More concretely, in the case of Qg x A4, the existence 
of the non-trivial centre implies that neutrino masses are generated via (£H£H)(ft 2n for some 
n > and the little group of ((ft) 2n is enlarged by the centre to (S, Z(Qs x A4)). A similar 
reasoning applies to the other candidate groups. Note that, so far we only investigated one 
flavon (ft in an irreducible representation, which does not apply in the more general discussion 
with multiple flavons (fti (or equivalently a reducible representation (ft). In this more general 
setup, the relevant combination of flavon VEVs contributing to the neutrino mass matrix can 
break the invariance again. Ultimately, the minimisation of the flavon potential decides which 
VEV alignment is achieved. 



3.2 General Group Extensions 

Let us have a closer look at the construction in the last section. In order to obtain the same 
flavour structure within G as within H, we demanded the existence of representations pi, 
which are directly related to the representations pf 1 of H. The representations p\ can be 
explicitly constructed using the surjective homomorphism from G to H, which we will denote 
by f : G -> H: 



14 The centre of a group, Z(G), is the set of elements, which commute with all elements of the group G, 
i.e. Z(G) = {x £ G \ gx — xg Vg £ G}. It forms a normal subgroup of G, i.e. Z(G) < G. 

15 The role of the unbroken subgroups in neutrino mixing has also been discussed from a bottom-up perspec- 
tive in [13]. 

16 The little group G/m is the subgroup of G leaving a VEV (</>} invariant, i.e. G/m = {g € G\g (<f>) = Wi- 
lt is also denoted by stabilizer subgroup or isotropy group. 
17 For model building of this type, see [26]. 
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Quotient Group H 


Order of G 


GAP 


Structure Description 


A A 


96 


201 


Z 2 .{Zl x A 4 ) 


144 


127 


Z 2 .(At x 5 3 ) 


192 


1017 


Z 2 .{D 8 x AO 


St 


96 


67, 192 


Zt.St 


144 


121, 122 


Z%.St 


192 


187, 963 


Z&.St 


192 


987, 988 


Z 2 .{{Zl x At) x Z 2 ) 


192 


1483,1484 


Z 2 .{Zl x Si) 


192 


1492 


Z 2 .((Zl x Z 3 ) x Z 2 ) 


r 


192 


1007 


Zl{Z% x At) 



Table 3: Candidate groups G up to order 200, which can not be written as semidirect product. 
The expression of the form N.H in the last column is the GAP notation of a central extension, 
i.e. N is a normal subgroup of G, which is contained in the centre of G, and H is the quotient 
group G/N H. Here, we explicitly choose N = Z[G) and therefore N.H = Z{G).G/Z{G). 
The candidate groups of order 200-500 can be found in Tab. 8. 

Hence, as soon as there is a surjective homomorphism £ : G —> H, there are representations 
Pi with the desired property. Therefore, is it enough to look for groups G and a surjective 
homomorphism £ : G —> H. This automatically implies the existence of a normal subgroup 
N = ker£ and a quotient group G/N = H. Thus, we are only dropping the condition that 
H is a subgroup of G. Actually, this type of extension is a general problem in group theory, 
which aims to find all possible groups G given two groups N and H, such that G/N = H. In 
the mathematical literature, this is denoted by short exact sequence. One example of such an 
extension is T' . At is not a subgroup of T', but At = T' /Z 2 . In T' models [14], the flavour 
structure of the lepton sector is essentially described by the quotient group T' /Z 2 = At 
and the additional group structure, i.e. the two dimensional representations 2j, are used to 
describe the quark sector. Hence, group extensions of the kind we described are not limited 
to the VEV alignment, but can be used more generally to lift properties of one group H 
to a larger group G, which addresses additional questions in flavour physics. Therefore, we 
propose to use these kind of constructions more systematically. 

In this article, however, we are mainly interested in a solution to the vacuum alignment 
problem, and therefore, we do not consider these other possibilities further, but perform 
another scan looking for groups solving the vacuum alignment problem and we relaxed the 
first condition of the previous scan to 

1. G/N ^ H with H being one of the groups At, T 7 , V 18 , 5 4 , A(27), 

while keeping the other conditions. It turns out that there are only candidates for At, T' and 
S4 up to order 1000. We collect all candidates up to order 200, which are not contained in 

18 We included T' in this scan, although T' is an extension of A4 via T' /Z2 — A4. However, the second 
condition excludes several candidates for T' , because the Z2 in T' /Z% = A4 is a subgroup of the N in the 
second condition. 
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the previous search for semidirect product groups, in Tab. 3 and present the candidates of 
order 200 - 500 in Tab. 8. 



4 Smallest Group: x 

In this section, we discuss a model of lepton masses and mixings based on the smallest 
semidirect product group Qg x A4 in the catalogue obtained in the preceding section. After a 
brief description of the group, we discuss why it is necessary to employ more than one faithful 
representation of the full group and build a model with this particle content. We then show 
that the most general scalar potential has the desired accidental symmetry and thus allows 
for the correct vacuum alignment. 

4.1 Group Theory 

While the A4 subgroup is presented by 

(S, T\S 2 = T 3 = (ST) 3 = 1) , (10) 

the quaternionic subgroup Qg (also known as D±, the double group of the dihedral group of 
order 4) is defined by 

(X,Y|X 4 = 1, X 2 = Y 2 , Y~ l XY = X- 1 ) . (11) 

The semidirect product Qg x A^ we are considering here is defined by the additional relations 
between the generators of Qg (X, Y) and A4 (S, T) 

SXS~ 1 =X, SYS~ l = Y- 1 , TXT~ 1 = YX, TYT' 1 = X . (12) 

An explicit matrix representation of these generators for the relevant representations is 
given in Table 4 and the character table is presented in Table 5. The Kronecker products 

3i x3 i = l 1 + l 2 + l 3 + 3 is + 3 iA (13a) 
5 

3;x3 j= ^3 k (i^j) (13b) 

3_i x 4j = 4 X + 4 2 + 4 3 (13c) 
4l x 4i = Ii 5 + 3 lA + 3 2s + 3 35 + 3 45 + 3 5a (13d) 
4i x 4 2 = l 2s + 3 lA + 3 2s + 3 35 + 3 45 + 3 5a (13e) 

show that if one uses the unfaithful triplet x ~ 3i to break A^ in the charged lepton sector 
and the four dimensional faithful representation (j) ~ 4_i in the neutrino sector, there are no 
dangerous cross-coupling terms of the form (4>4>)^_^ (xx) 3 ^ e ^ c - allowed by the symmetry that 
would forbid the required VEV alignment. 
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v° 




-1 
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Table 4: Relevant Representations of Qs x A4 in some basis. The first 4 representations are 
the unfaithful A4 representations the leptons are assigned to (therefore X = Y = 1). The 
last representation is used to break A4 in the neutrino sector. Note that this representation 
is double valued, i.e. X 2 = Y 2 = — 1. Here uj = e l27r / 3 . 



The relevant operators for the generation of the lepton masses are £xf c H with / being e, 
li or r as well as £HIH and iHlHc/} 4 for neutrino masses 19 . 

Unfortunately, the most general VEV configurations 4> ~ (a> b, —b) that break the group 
to the Z2 subgroup generated by S cannot be realised in the flavon potential 20 



due to the relation 







dV A 



l x +ai 



i=2,3 



3: 



b 2 or a6 



a6(a 2 — b 2 )(ot2 + 03) 



(14) 



(15) 



lead to a restoration of symme- 



The achievable VEV configurations with o 2 

try in the operator {U)^_^ (^ 4 )3^ ^ na ^ generates the {U)%^ entry in the mass matrix and 



\21 



consequently it vanishes in the vacuum, ((<^ 4 )3^) ~ ab(a 2 — b 2 )' 

This type of model is also not so interesting from a general point of view, as it shares 
a couple of unpleasant features with the model of Babu and Gabriel [21] when viewed as an 
effective field theory: 

• the off-diagonal entries in the neutrino mass matrix, generated by (IHlHcf) 4 ), would be 
of very different order than the diagonal ones generated by the operator (IHIH). To 
satisfy neutrino data, the two entries have to be of almost the same size, though. 

-x, r 



19 Here we again assume the discrete Z2 symmetry \ 
neutral fermion sector. 



20 The operator 
of the other operators. 



— f c to separate the charged from the 



3 ' (^0)3 i which one would naively expect, can be expressed as a linear combination 



21, 



If one introduces a soft-breaking term that conserves the Z2 subgroup generated by S, Vs = a (4>i(f>2 + $3 4>4) 
in the potential, the minimum with a 7^ b can then be realised. We do not pursue this option further here, as 
we are interested in genuine spontaneous symmetry breaking. 
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1 T SYX SY Y 2 T 2 TY S SX X STYT 



It 

—1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


—2 


1 


Co' 


1 


1 


1 


Co- 2 


Co' 


1 


1 


1 


uj 2 


la 


1 


u 2 


1 


1 


1 


UJ 


Co' 2 


1 


1 


1 


UJ 


3l 


3 




-1 


-1 


3 






-1 


-1 


3 




3 2 


3 




3 


-1 


3 






-1 


-1 


-1 




3 3 


3 




-1 


3 


3 






-1 


-1 


-1 




3 4 


3 




-1 


-1 


3 






3 


-1 


-1 




3 5 


3 




-1 


-1 


3 






-1 


3 


-1 




4l 


4 


1 






-4 


1 


-1 








-1 


4.2 


4 


a; 2 






-4 


Co' 


-cu 2 








-UJ 


4 3 


4 


Co' 






-4 


UJ 2 


-Co' 








-UJ 2 



Table 5: Character table of Qs x A4. The first line are representatives of the different 
conjugacy classes. Zeroes in the character table are denoted by a dot . and uj is the third root 
of unity uj = e 27r1 / 3 . 

• as (£H£Hx 2 ) is allowed and of smaller dimension than (IHlHcj) 41 ), tri-bimaximal mixing 
is not a leading-order prediction of the model. 

All of these issues can of course be cured by introducing a UV completion that does not 
confirm the effective field theory prejudices. 

However, we restrict ourselves to natural solutions within effective field theory. To solve 
all of these problems, we will discuss a model with two flavons in the neutrino sector, </>i ~ 4^ 
and 4>2 ~ 4j_, where an additional symmetry forbids the allowed term x ■ (4'i4'2)s-^ that could 
disturb the VEV alignment between the two sectors. We identify this symmetry with the 
one that separates the charged lepton from the neutral lepton sector, i.e. we postulate the 
additional Z4 symmetry I — > i£, f c — > —if c and 4>2 -4 —(p2, where / denotes e, [x and r. One 
can think of this Z4 symmetry as a discrete version of lepton number with being doubly 
charged under this lepton number. 

4.2 Model and Lepton Masses 

Finally, we present a model based on the symmetry group Qs x A4 augmented by the auxil- 
iary symmetry Z4 introduced at the end of the last section. The leptonic and scalar particle 
content is given in Tab. 6. As advertised, for the standard model leptons, we use the un- 
faithful representations 1^ 2 3 an( ^ — \, that transform as irreducible representations under 
the subgroup A$. In the charged sector we use the unfaithful representation \ ~ 3^ and 
the charged lepton sector is thus analogous to the usual construction in an model. In the 
neutrino sector, we introduce the real flavons ~ 4^. 

To keep the discussion simple, we use an effective field theory description. To lowest order, 
the charged lepton masses arise from the operators 

4 5) = Veiixh^H/A + y^x)l 3 ^ C H/A + y^l^H /A + h.c. , (16) 
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with H = i<T2-ff*, and the neutrino masses are generated from the effective interactions 

4 ? ) = x a {tmH) Xi {<h<h)i x /A 3 + x^miH)^ ■ (0i0 2 )3 l /A 3 + h.c. . (17) 

The notation should be self-explanatory and the relevant Kronecker products are given in 
appendix A. We will show in the next section that the vacuum configuration 



(x) = (v',v',v') 



t\T 



1 



(a, a, b, —b) 1 



1 



(c, c, d, —d) J 



(18) 



with v' ,a,b,c,d G M, can be obtained as the global minimum of the most general scalar 
potential. This configuration gives ((^1^2)3^) = \{bc— ad, 0, 0) T and ((^1^2)1^) = \{ac+bd) 
and it is of course non-unique as there are many more physically identical patterns that can 
be obtained by acting with the group generators on this vacuum. This VEV configuration 
breaks the flavour symmetry to the Z 2 subgroup generated by S. There are also physically 
inequivalent minima of the potential that break to the Z2 subgroups generated by SY and 
SYX which lead to the same structure ( (^1*^2)3-) oc (1,0, 0) T . We will comment on this 
more in sec. 4.3. 

The leading-order mass matrices are given by 



M E 



v'v 



V2K 



' Ve 
\ Ve 



with 



x a -^(ac + bd), 



Vt 



x d^(bc-ad) 




(19) 



and 



(H) 




(20) 



The mass matrices can be diagonalized by 



Un 



1 

71 



fi 1 
1 
1 



2 




(21) 



such that U^M, 



-^diag(y e ,y^,y T ) 22 , UjM u U u = diag(a + d, d, -a + d) and the 

resulting mixing matrix Umns = UqU u is given by the HPS matrix (1). This construction is 
of course completely analogous to the usual A4 models[10] and it is well-known that the mod- 
erate tuning \d\ ~ \d\ is needed in order to accommodate the correct neutrino spectrum[27]. 
However, in the usual models the contributions to a and d stem from completely different 
VEVs, while in our model both stem from VEVs of the same fields, and a similar order of 
magnitude might therefore be considered more natural. Indeed, in the numerical minimisation 
of the potential, we found a tendency for a similar size of the two <f> contractions. 



22 The charged lepton mass hierarchy can be explained by a Froggatt-Nielsen U(l) symmetry in the usual 
way. 
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particle 


ST/(3) C 


SU{2) L 


U(l) Y 


Q 8 x At 


Z 4 


t 


1 


2 


-1/2 


3l 


i 


e c + [i c + r c 


1 


1 


1 


I1+I2+I3 


— i 


H 


1 


2 


1/2 


ll 


1 


X 


1 


1 





3l 


1 


01 


1 


1 





4l 


1 


02 


1 


1 





4l 


-1 



Table 6: Particle content of the minimal model with the correct spontaneous symmetry 
breaking. All fermions are left-handed Weyl fermions. 



Let us also comment on the fact that we have to employ two real copies of the faithful 
representation 4^. This exactly matches the numbers of degrees of freedom of one complex 
A4 triplet and one complex singlet, which is commonly used (e.g. [16]). Here, we do not have 
to introduce additional degrees of freedom to obtain the correct vacuum alignment and we 
thus think it is an attractive and economical model. 

The effects of higher order operators can be found in sec. 5. 

4.3 Scalar Potential 

Here, we demonstrate that the pattern of vacuum expectation values we used in the last 
section can be obtained as the minimum of the scalar potential. The most general scalar 
potential invariant under the flavour symmetry is given by 

F(x,0i,02) = V x (x) + ^(0 1 ,0 2 ) + y mix ( X ,0i,02) (22) 

with 

K?i(0i,02) =^1(0101)1-, +ai(0i0i)| + ^2 °*(0i0i)3. " (01 0i)3_. 

i=2,3 1 1 

+/i!(0202)l 1 +/3l(0202)| + ^2 A(0202)3. ' (0202)3- 

i=2,3 1 1 

+7i(0i0i)l 1 (0202)l 1 + ^2 7i(0i0i)3. • (0202)3. 

i=2,3,4 1 1 

v x (x) = ^l{xx)\ 1 + pi{xxx)\ 1 + AiCxx)^ + a 2 (xx)i 2 (xx)i 3 

Vmix(X, 01, 02) = (13(0101)^ (XX)^ + (23 (02 02^ (XX)^ (23) 

Note that, by construction, there are no non-trivial couplings between the x an d breaking 
sectors that would disturb the vacuum alignment. The potential thus has an 'accidental' 
HQs x A4) x A4] x Z4 symmetry. This symmetry is explicitly broken to (Q$ x A4) x Z4 
by the couplings to leptons and by higher dimensional operators in the potential. As the 
accidental symmetry is discrete, there is no pseudo-Goldstone boson, as can easily happen in 
constructions of this type [21]. 
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Let us now demonstrate that this model does not suffer from a vacuum alignment problem. 
At first, we discuss the possible minima of the potential focussing on the little group in the 
neutrino sector, i.e. the subgroup, which leaves the VEV invariant. If there is a minimum, in 
which the symmetry generator Q £ G is left unbroken, i.e. Q (fap) = (^1,2) > there obviously 
are degenerate minima (^fa t 2^ = g (fa,2) that leave gQg -1 unbroken, with g £ G. The 
physically distinct minima are therefore characterised by the conjugacy class (es) G ■ Qi = 
{gQig~ 1 \g £ G} of the group element(s) Qi. Obviously, only conjugacy classes with an 
eigenvalue +1 can lead to a non-trivial little group. For the four dimensional representation 
4 1 , there are five such classes which are represented by 1, S, SY, SYX, T as well as T 2 . The 
groups generated by T and T 2 are identical. For the three dimensional representation 3^, 
where X and Y are represented trivially, all conjugacy classes have an eigenvalue +1 and can 
lead to a non-trivial little group. The relevant little group in the neutrino sector is the one 



of {(fafa)^)- 

In the following, we will firstly discuss the possible little groups of (fa) and then its 



implications for the little group of ((fafa)^ )■ There are three physically distinct minima 



of fa, that preserve a Z 2 subgroup: 

• {fa) = ^(«, a, b, —b) T results in the little group {S), 

• {fa) = (0,a,6,0) T in (SY) and 

• = ^(.~a,b,-a,b) T in (SYX) . 

In addition, there is one preserving a Z3 subgroup: 

• (fa) = A=( a ,a,a,b) T preserves (T) (as well as <T 2 } = (T)). 

Obviously, there are also minima leading to little groups, which are generated by more than 
one generator. For example (fa) oc (1, 1, 1, — 1) T preserves (S, T) = A4. The same discussion 
applies to fa. The little group of ((fafa) 3l ) contains the intersection of the little groups of 
(fa) and (fa). 

In the following, we will concentrate on the three little groups (S), (SY) and (SYX), 
which we listed above. If both (fa) and (fa) preserve the same Z 2 subgroup, we obtain 
{{fafa)^) = \(bc - ad,0,0) T and {(fafa) ll ) = \(ac + bd) due to X = Y = 1 for the 
SM representations (lj and 3^) with a, b being the VEVs of fa and c,d, the corresponding 
ones of fa. Thus, it is impossible to distinguish these minima from low-energy neutrino 
phenomenology at the leading order. They are, however, physically distinct, since (fafa) as 
well as (fafa) are different and they lead to different mass spectra in the scalar sector. The 
minimisation conditions for the VEVs are: 



a (a+ (a 2 + b 2 ) + a_ (a 2 - b 2 ) + 7+ (c 2 + d 2 ) + 7 _ (c 2 - d 2 ) + Ui) + Tbcd = (24a) 

b (a+ (a 2 + b 2 ) - a_ (a 2 - b 2 ) + 7+ (c 2 + d 2 ) - 7 _ (c 2 -d 2 )+Ui)+ Tacd = (24b) 

c (0+ (c 2 + d 2 ) + /?_ (c 2 - d 2 ) + 7+ (a 2 + b 2 ) + 7 _ (a 2 - b 2 ) + U 2 ) + Tabd = (24c) 

d (P+ (c 2 + d 2 ) - f3- (c 2 - d 2 ) + 7+ (a 2 + b 2 ) - 7 _ (a 2 - b 2 ) + U 2 ) + Tabc = (24d) 





(24e) 
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where the equations have been rescaled to eliminate constant overall factors and with the 
shorthand notations 

Ui = \(£ + V3( l3 v' 2 ) for * = 1,2, U 3 = 2nl + ( 13 (a 2 + b 2 ) + ( 23 (c 2 + d 2 ) 

and 

r u = ^, t- = ^ for z = a 'P 

\ 7+ = ^, 7 _ = ^ and T=% 

_ ^71+73 _ _73_ J p = 72+74 

e+ = S 1 iV% ) e- = ^ for £ = a,0, 

_ V371+72 _ 73 +74 j p _ 72 

The first four equations result from the derivatives taken with respect to the components of 
4>i and cp2- The eleven minimization conditions, corresponding to the 11 real scalar degrees of 
freedom, thus reduce to just five equations for five unknowns and there generally is a solution. 
Thus, there is no vacuum alignment problem in this model. We have checked numerically 
that this is the global minimum for a region of parameter space. 

Note that the equations for v' and a,b,c,d essentially decouple and the contribution to 
the other one, can be reabsorbed in the mass term. They are invariant under symmetries 
(a, c) ++ (6, d), (a, b) — > —(a, b), (c, d) — > — (c, d) as well as (a, b, cti, U\) -f+ (c, d, C/jj), which 
are inherited from the symmetries of the potential. 

There are also minima breaking all symmetries. Generally, for each minimum, there are 
95 additional minima with the same value y| m in> which are connected by a group transforma- 
tion. This multitude of minima makes an analytic treatment unfeasible. Therefore, we have 
performed a numerical study. We have varied all parameters in the range [—4,4] and only 
found minima corresponding to the ones given above. In particular, we have not found any 
minima that leave a larger symmetry group intact (except for the minimum with vanishing 
VEVs, which does not break the group). Each of these minima can be realised as global mini- 
mum of the potential, which we checked in the random number scan we performed. However, 
it was impracticable to determine the parameter regions of each global minimum. 



5 Higher Order Corrections 

The results presented above are corrected by higher order operators. Here we discuss the 
next-to- leading order corrections. Let us briefly comment on the magnitude of the scale A 
under the assumption that all operators are suppressed by the same scale. 23 If we require a 

23 0f course, this assumption does not have to be true for e.g. a UV completion where the charged lepton 
mass operators are generated by vector-like fermions and the neutrino mass operators are generated by a 
see-saw. 
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perturbative value for the r Yukawa coupling, y T < 4tt, this translates into [16] 



v' 

— > 0.002. 
A 

Furthermore taking m u ~ 0.05 eV and assuming couplings of order one, x a ^ ~ 0(1), we find 

A PS 6- 10 14 GeV with u = a,b,c,d 

and the natural cutoff values are therefore 2 • 10 9 GeV < A < 6 • 10 14 GeV assuming all VEVs 
to be of a similar size v' ~ a ~ b ~ . . . , but it can easily be in the TeV region for moderately 
small couplings in the UV completion. 

5.1 Corrections to the charged lepton mass matrix 

The next-to-leading order correction to the charged lepton mass matrix takes the form: 

4 6) = y'Mxxh^i^H/k 2 + y'.Wxxh^l^H/A 2 + y'Mxxh^i^H/k 2 + h.c. . 

(25) 

As these operators can be obtained by replacing \ by (xx)s-^ m ^q- 0-®) an< ^ 

((xx)^) = v '(x), 

they do not introduce a new structure in the charged lepton mass matrix [28], but merely 
renormalise the leading contribution. Note that there are no other contributions at this level, 
since 4>i^>i does not contain 3 \ by construction. Operators with new structures are suppressed 
by 1/A 3 . 

5.2 Corrections to the neutrino mass matrix 

The next-to-leading order operators contributing to the neutrino mass matrix are given by 

A 4 4 8 ) = . c (fflffl) l2 (^to)i 3 +x 6 (£M)i 3 ((M 2 x)i 2 +* h (*ff£ff) il fai6 2 x)l 1 + 

+ (eHM)3 i ■ [xexfal&hj +xf(x- (fafafo^S + Xg(X' (fafa)^ )a] + h.c. , (26) 

where (...)s denotes the symmetric contraction and (...)a the antisymmetric one. These 
operators perturb the mixing matrix and their effect will be discussed in section 5.4. 

5.3 Corrections to the Scalar Potential 

Corrections to the potential arise at dimension five: 

2 4 fi(LM) 



V(5) = E E • (0l0l)3. • (Mm)3. + 

L,M=1 i,j=2 1 J ' % 



+ X + 4 3) (^1^1)1! + ^(^hj (27) 
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where all parameters are real and <5- • 



(LM) 



for i > j. Upon minimisation, these interactions 



lead to a shift in the vacuum expectation values of the form: 



(x) = (v + 5v[,v' + 5v' 2 ,v' + 5v' 2 ) T , 



1 
1 

^2 



(a + Sai, a + 5a 2 , b + 5a 3 , —b + 5a^) T , 
(c + 5b 1 ,c + 5b 2 ,d + 5b 3 , -d + 56 4 ) T 



Generically, the magnitude of these shifts will be suppressed by one power of A, 

Su u 
u A' 



(28a) 
(28b) 

(28c) 
(29) 



where u denotes a generic vacuum expectation value. The VEVs of X2 and X3 stay equal at 
next-to-leading order, i.e. (X2) — (X3) = C(l/A 2 ) and (^3) ~ (X2) = 8v' 2 . To calculate the 
correction to neutrino masses, the following shorthand notations for the shifts in the vacuum 
expectation values are useful: 



and 



f a (Sbi — 8b 3 ) + c (5a 3 — 60,4) — d (Sat + 5a 2 ) + b (5bi + 5b 2 ) 
a (5b 3 + Sbi) — c (Sas + £04) + d (Sai — 5a 2 ) + b (Sb 2 — Sb\) 
y a (Sb\ — Sb 2 ) + c (Sa 2 — Sai) — d (Sa 3 + £04) + b (Sb 3 + Sb^j 



( 5<&i 

<5$ 3 

(30) 



(a (6h + 5b 2 ) + c (Sai + Sa 2 ) + d (Sa 3 - 5a 4 ) + b (Sb 3 - Sb A )) = <5$ • (31) 



To get a feeling for the size of the deviations from the leading order vacuum alignment, 
we have performed a numerical minimisation of the potential for a number of random values 



for the potential parameters. We found it instructive to plot ^S-tT' against 



A ' 



where m 



denotes any of the leading-order VEVs and Sui any of the deviations. Fig. 1 shows the VEV 
deviation scales plotted against the ratio u/A. The corrections are small for small u/A. 

5.4 Corrections of Masses and Mixings 

To next-to-leading order, the charged lepton matrix Me is modified from Eq. (19) by 



/ Sv[ 



5M E 



AV2 



5v'o 



Un 



\ 







Sv'o 



( Ve 

y M j + 
V y r 



vv 



A 2 v/2 





V y' T 











(32) 



In the neutrino sector there are also new structures. The corrections to the neutrino mass 
matrix can be parametrised as 
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Figure 1: Vacuum shifts maxjjttl'l mcm ced by higher dimensional operators as a function of 
max^|Mfc| ran( j om ^y c hosen potential parameters of order unity. All points correspond to 
phenomenologically viable data points. 
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As the leptons only transform under the A4 subgroup of the model, the neutrino phe- 
nomenology runs exactly parallel to the A4 case. The effects of the operators a,...,f have 



been studied in [29] where it has been shown that a sizeable deviation from sin 2 



'13 



is pos- 



sible without introducing large corrections to the other mixing angles. Recently it has been 
shown that sin 2 #13 > 0.1 is possible for c/d > 0.25 in the case of normal mass ordering [18]. 

We performed a scatter plot in order to get an idea of the size of the corrections from 
higher dimensional operators. For a collection of tree- level parameters of order unity, we have 
varied the higher dimensional parameters (27) of the potential in the range [0.5, 1.5] and the 
dimensionless parameters in the corrections to the lepton masses in Eq. (26) have been taken 
to be of the same order as the leading order contributions. The suppression scale A has been 
varied in a wide range. 

In Fig. 2, the resulting scatter plots are shown, where all data points lie within the 3<r 
limits of the global fits cited in the introduction. As can be seen from Fig. 2(d), for u/A > 0.05 
there are points that deviate from tribimaximal mixing in the right way to be compatible with 
the recent measurements from T2K. 
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Figure 2: Scatterplot of Mixing Angles. To illustrate the typical size of corrections to the 
mixing angles, we have performed a scatter plot. We took all dimensionless scalar potential 
couplings to be of order one and varied the ratio of the mass parameters in the potential such 
that the ratio of the VEVs and cutoff- scale is smaller than one. All dimensionless parameters 
that modify the neutrino and charged lepton mass matrices are taken of the same order as 
the leading order parameters. All points lie within the 3 a range mass of the mass and mixing 
parameters. The mixing angle sin 2 9\2 is varied more than the the other two mixing angles. 
We have used the MixingParameterTools [30] package to extract the mixing angles. 



Allowing for couplings considerably smaller than order one in V^ 24 , the VEV corrections 
become dominant and e and / are the main corrections to the neutrino mass matrix. 
This is shown in Figure 3 and is in agreement with the result [29]. Note that the values are 
roughly along a diagonal line, i.e. e and / are similar in size, but have a different relative 
sign. 

5.5 Cosmological Implications of Accidental Symmetries 

Let us briefly comment on possible cosmological implications of the unbroken remnant Z<i and 
Z3 symmetries in the neutrino and charged lepton sectors, respectively. Due to the accidental 

24 For details, please consult the Mathematica notebook published as a supplement together with the Math- 
ematica package Discrete described in sec. 8 
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Figure 3: The matrix entries e and / are the dominant corrections to the neutrino mixing 
matrix in the case where the corrections from the VEVs of <j)\ and <p2 dominate. Here we 
show the correlation between the two quantities. The blue, violet, yellow and green points 
correspond to values of sin 0? 3 in the ranges [0,0.005], [0.005,0.01], [0.01,0.02], and [0.02,0.03], 
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A4 symmetry in the scalar potential these symmetries are accidental symmetries of the theory 
only broken by higher dimensional operators. 

Let us discuss the situation where J is the lightest scalar odd under the unbroken Z2 
symmetry generated by S, e.g. J = ^ ((^1)3 + (^l^)- It can then decay into neutrinos 
through the effective interaction 

£ = -T^gjv^JviVj + h-c- (35) 

with a lifetime roughly given by 

167T u 2 o / u \ ( u \ , . 

r( ^^~^~ 4 ' rf H^)(«s) (36) 

for m u = 0.05 eV and u being a generic flavon VEV . Depending on the model parameters, 
this decay time can be problematic. If the lifetime is larger than the age of the Universe, 
J becomes a dark matter candidate. A large lifetime naturally occurs, if J is a pseudo- 
Goldstone boson [31], which leads to mj/u <^ 1. Pseudo-Goldstone bosons often appear in 
these constructions. For example, the tree-level scalar potential of the next-larger group in 
Tab. 2, T' xi A4, has the large continuous accidental symmetry Sp(4). 

However, in general, there is also the decay channel via higher dimensional operators in 
the scalar potential, which couple J to the (£) -breaking VEV of \i e -g- by operators of the 
type 4>f ■ xH^H. It will generically be the dominant decay process in the model outlined above 
and result in much shorter lifetimes of 

T ~ 16^ -3.3. 10- (">) ('^Y 7 ( 10l2 ?° V ) . 
u 8 V u J V0.01 J \ A J 
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ensuring that any potential abundance of J will decay before big bang nucleosynthesis. In the 
model by Babu and Gabriel [21], these higher dimensional operators are absent and therefore 
this decay through neutrinos is the only decay channel, which poses a potential problem for 
these models. 



6 See-saw UV completion 

The neutrino sector of the effective theory outlined above, may be UV completed by intro- 
ducing the left-handed Weyl spinors N, S2 and S3 that transform under (Qg x A4) x Z4 as 
N ~ (3]_, — i), 5*2 ~ (42, i) and S3 ~ (43, — i). S2 and S3 can be combined in a Dirac spinor. 
This leads to the following new interactions in the Lagrangian 

C = X£ N £HN + xn2NS24>i + XJV3-/VS302 + mS2S3 + X23S2S3X + h.c. , (37) 

where the contraction of each operator is uniquely determined by the group theory of Q$ x A4. 
The neutral fermion mass matrix is then schematically given by 

/ x m (H) \ 

1 ... x N2 {(pi) x N3 ((p 2 ) 

2 m + x 2 3{x) 

\ / 

in the basis (u, N, S2, S3). In the following, we assume that the direct mass term is larger 
than the mass terms generated by VEVs. Therefore, we are in the seesaw regime, which has 
been firstly studied for gauge singlets in [32] and in more generality in [33]. Hence, the masses 
of the singlets N are generated 

/ A 

A B ] with A = -2(ac + bd) and B = iV3(bc - ad) . (39) 



(38) 



XN2XN3 

m N = 



m 








This particular form has been denoted linear see-saw [34]. The light neutrino masses are 
generated via a standard see-saw [32]. Hence, the operator X23S2S3X does only enter at next- 
to leading order. Alternatively, it is possible to forbid it together with all next-to leading order 
corrections, which have been discussed in the previous section, by introducing an additional 
Z2 symmetry \ ~ * ~ X anci f c ~^ ~f° with / being e, fi or t. The neutrino mass matrix is 
then given by 

m v = xf A rU 2 m^ v 1 (40) 

This can also be seen from Figure 4. This matrix is diagonalized by U v : U^m v U v = 
diag( g ^ , B ^_ A )- However, there are two degenerate eigenvalues as the relative phase 
of A and B is given by tt/2. This can be solved by adding another copy of S2 or S3, for 
example, lifting the degeneracy. 

The charged lepton mass operators can be generated in the same way as in [21] by intro- 
ducing additional states that have masses allowed by EW symmetry and mix with the SM 
states after EW symmetry breaking. 
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Figure 4: Neutrino masses in the UV completion. 
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Table 7: Chiral Superneld Particle Content. 



7 Supersymmetrization 

Supersymmetrization of the model is rather straightforward. One only has to ensure that 
there are no flat directions in the potential: 

V = V SVSY + y soft (41) 

with 

VSUSY = / J 

i 

where W denotes the superpotential and (f{ is any of the fields in Table 7. V so ft contains all 
supersymmetry-breaking soft terms invariant under the flavour symmetry. 

As there is no cubic invariant containing the fa 2 fields only and the quadratic term 
2 = Yli 01 2 i s SO(4) 2 invariant under the individual rotations of fa and fa, we have to 
add the singlet S and the triplet x ~ 3_2, to get a superpotential without flat directions in 
the cubic terms and without a continuous accidental symmetry. We thus have the schematic 
superpotential 

w = s{4>\ + 4 + x 2 + x 2 )^ +s 3 + s 2 + s + 4> 2 1 + ( f> 2 2 + x 2 + x ;i + x(4>\ + 0i)3 2 + x 2 + x 3 - 

(43) 

We have studied the potential resulting from this superpotential and the most general soft- 
breaking terms and we have found a portion of parameter space with the right vacuum 



dW 



dcpi 



(42) 
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alignment, with non- vanishing VEVs for both the singlet and triplet contractions of the 
product 0i02- The neutrino mass operators are again given by 

W D XaiWeH^ifafch^A 3 + x d (^^if u )3 i (0 1 2 )3 i /A 3 . (44) 

As in the non-SUSY model before, the the on-and off-diagonal terms of the neutrino mass 
matrix, which have to be quite close to each other in magnitude, are generated by VEVs 
of the same fields. The additional scalar field \ couples to leptons only on next-to next-to 
leading order and it is thus not problematic. 

Details can be found in the Mathematica notebook accompanying this paper, which can 
be downloaded from the webpage of the Mathematica package Discrete. We do not give 
the details here, as it is somewhat out of the main focus of the paper, but we have checked 
that there exist parameter values for which the global minimum of the potential has the 
correct vacuum alignment for the most general softly broken supersymmetric potential. The 
symmetry breaking is also complete, i.e. there are no flat directions left as is the case in the 
model by Altarelli and Feruglio [28], where one introduces driving fields with Ur(1) charges of 
2. It has been shown that the inclusion of soft-breaking terms in this model is problematic, as 
it generically leads to flavour violating VEVs of auxiliary fields [35] , unless there is a solution 
to the SUSY flavour problem in terms of gauge mediation [36] or via a separate mechanism 
(see e.g. [37]). 

It would be interesting to come back to the SUSY VEV alignment problem and search for 
groups that do not need additional scalar fields to break accidental symmetries. 

8 Discrete — Mathematica Package 

Discrete is a Mathematica package with several useful model building tools to work with 
discrete symmetries. The main features are 

• the calculation of arbitrary Kronecker products, 

• an interface to the group catalogues within GAP [22], e.g. the SmallGroups [ ] library 
with all discrete groups up to order 2000 (with the exception of groups of order 1024) 
and many more. 

• calculation of Clebsch-Gordan coefficients. They are calculated on demand and are 
stored internally, in order to improve the performance. 

• the possibility to reduce covariants to a smaller set of independent covariants. 

• the documentation is integrated in the documentation centre of Mathematica. 

Discrete can be downloaded from http : //projects .hepforge . org/discrete/. It has been 
tested with Mathematica 8 running on Linux as well as MacOS, but we expect it to run on 
older versions of Mathematica as well. 
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It requires a working installation of GAP [ ] as well as the GAP package REPSN [38]. GAP in- 
cluding all its packages can be downloaded from http : //www . gap-system . org/. On Debian- 
based Linux-distributions, it can be directly installed via the package management. 

In the following, we are presenting a short example of the abilities of Discrete and refer 
the interested reader to the documentation and the example notebook within the package. 
For simplicity, we are choosing Aa and only calculate the renormalizable part of the flavon 
potential of <fi ~ 3. For brevity, we have shortened the output. The omissions are denoted by 
dots. 

Needs ["DiscreteModelBuildingTools"]; 
A4=MBloadG APGroup [" AlternatingGroup(4)" ] ; 

Dimensions of irreps: 
12 3 4 
1113 

phi=MBgetRep Vector [ A4,4,p] 

$ Assumptions= Variables [phi] € Reals; 

{{},{MM{pi,p2,p3}}} 

V=MBgetFlavonPotential[A4,phi,4,h] 

h2nl (pi 2 + p2 2 + p3 2 ) 1 , „ , 9 . 2 

h3nl pi p2 p3 + K - ^ + -h4n2 (pi 2 + p2 2 + p3 2 ) 

y3 3 

1 , 1 a 99 a 9/9 9\ \ h4n3 (p2 2 p3 2 + pi 2 (p2 2 + p3 2 ) ) 

+ -h4nl (pi 4 + P 2 4 - P 2 2 p3 2 + P 3 4 - pi 2 (p2 2 + p3 2 )) + ^-^ ^-^ V —^- 

3 y3 

After loading the package in line 1 and loading the group Aa from GAP, we define a field 
phi in the third line in boldface, which transforms as triplet of Aa. The 4 denotes the triplet 
in the list of representations and the last argument determines how the components of phi are 
denoted. Furthermore, we declare all components of phi to be real. MBgetFlavonPotential 
returns the flavon potential of phi up to fourth order as specified in the third argument and 
the couplings start with h. The first number in the name of the coupling denotes the order 
and the second one enumerates the couplings of a given order. 

Part of the calculation for Q$ x Aa_ is included in Discrete as example. However, we 
recommend to start with the tutorial included in Discrete, which introduces and explains 
most functions. 

Recently, a Mathematica package has been presented that allows one to calculate the 
group invariants formed from the three dimensional representation for most finite subgroups 
of SU(3) with order smaller than 512 [39]. 



9 Conclusions 

In this paper, we have revisited the long-standing problem of vacuum alignment in models 
with a discrete flavour symmetry. In such a model, in order to obtain the correct pattern for 
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the mixing angles, it is generally necessary to break the flavour group in a specific way to two 
different subgroups. This vacuum alignment, however, cannot be realized as a minimum of 
the scalar potential due to non-trivial couplings between the two sectors responsible for the 
breaking to the different subgroups. 

We have performed a systematic scan of all discrete groups with less than 1000 group 
elements. For each of the flavour groups H = A4, T7, S4, T' and A(27) of the SM fermions, 
we have identified a number of candidate groups G with G/N ~ H and a three dimensional 
representation x inherited from H. We further required the existence of a faithful representa- 
tion (f), such that there appears an accidental symmetry G x H in the renormalizable part of 
the scalar potential. The flavon x, responsible for symmetry breaking in the charged lepton 
sector, and the SM fields essentially only transform under the group H, thereby preserving 
the mixing angle predictions of H. The flavon <p breaks the symmetry in the neutrino sector 
and its product (j) x (j) does not contain any of the representations of H . This is a necessary 
condition for the accidental symmetry, as the term (j) 2 can - infamously - not be forbidden 
by an internal symmetry with unitary representations acting on (J). The additional symmetry 
thus forbids the dangerous cross-couplings, i.e. there is only the trivial coupling via the total 
singlet between <f> and x a t the renormalizable level. The accidental symmetry is then broken 
to G by couplings to fermions and other higher dimensional interactions. 

Having identified a list of possible groups, we built an explicit model using the smallest 
semidirect product of the candidates, Qg x A4, as flavour group. We have used two real scalar 
copies <pi and (j)2 of the faithful representation 4]_, where the triplet contraction (<^i</>2)3^ 
couples to neutrinos and thus plays the role of 4>s in the model of Altarelli and Feruglio [16]. 
The accidental symmetry is protected by an additional Z4 separating the charged and neutral 
lepton sectors. We have explicitly shown that the potential has the desired vacua and that it 
does not lead to unwanted pseudo-Goldstone bosons, i.e. the symmetry breaking is complete. 
We have further discussed the influence of next-to-leading order higher-dimensional operators 
on masses and mixings. 

As a direction of future work, it would be interesting to study a model where the flavons 
X and 0i 2 transform in the same way as the Higgs field under electroweak symmetry, which 
would move the flavour breaking scale to the electroweak scale and make it testable. We 
think that this model is quite well suited for this study as there are two Higgs fields in the 
Weinberg operator and only one in the Yukawa couplings, as fits nicely with the structure of 
our model. 

Let us conclude by a brief comparison with other schemes of obtaining the correct vacuum 
alignment. Counting degrees of freedom of the effective theory, our model has the same 
number of degrees of freedom as the minimal model without any mechanism for vacuum 
alignment [ ]. While we here do not have to add any degrees of freedom, the solution 
of the VEV alignment problem with an XJ(1)r symmetry as well as a brane constructions 
require a plethora of additional degrees of freedom in the form of driving fields or KK modes, 
respectively. When compared with the model of Babu and Gabriel based on the wreath 
product of S3 with A4 [21], our model not only has a substantially lower number of degrees 
of freedom, but it also works as an effective theory, as (£H£H)q and (£H£H)-^ are created 
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on the same order, which is not the case in their model. 

Finally, it might be worthwhile to look into other extensions of flavour groups used in the 
lepton sector to address for example the quark sector. One prominent existing example is 
the extension of to T' , which enables to describe the lepton and quark flavour structure 
simultaneously. We expect that our approach described in sec. 3 will be a useful tool for 
model building in this direction. 
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A Clebsch-Gordan Coefficients 

In this section, we present the Clebsch-Gordan coefficients, which are relevant for the discus- 
sion. 

A.l A 4 

The only non-trivial Kronecker product of A4 is given by 



where the indices S and A indicate whether the representation is in the symmetric or an- 
tisymmetric part, respectively. The corresponding Clebsch-Gordan coefficients, which have 
been computed using [40], are 



3 x 3 = l x + 1 2 + 1 3 + 3 5 + 3 A 



(45) 




1 



(46) 



where (ai,a 2 ,a 3 ), (h,b 2 



63) ~ 3. 



29 



(a) G/N ^ A A 



\G\ 


GAP 


240 


108 


288 


924 


336 


131 




618 


5809 


384 


5856 


18216 




20112 




432 


262 


480 


964 


1041 







(b) G/N 


s 4 




\G\ 


GAP 




X \J £j 


X WO 






288 


4DD 
o^± i 


844 
on 1 } 

yuo 


845 


846 


oou 


1 1 5 

X X (J 


X X \J 








582 


5614 


5705 


5708 




5713 


5714 


5728 


5733 


384 


18028 


18029 


18042 


18043 


18044 


18045 


18048 


18102 




18117 


18120 


18130 


18143 




20069 


20073 






432 


240 


241 






480 


257 


961 


967 


968 


969 


970 


1020 





(c) G/N^T' 



\G\ 


GAP 


288 


409 


384 


5845 


480 


266 



Table 8: Candidate groups G of order 201 — 500. |G| denotes the order of G. The groups up 
to order 200 are listed in Tab. 3. Details of the groups may be accessed using the computer 
algebra system GAP by using the command SmallGroup (Order , GAP) . 



A. 2 Q s X A 4 

The product of two triplets 3j x 3j is described by the same Clebsch-Gordan coefficients as the 
one in A±. They are shown in Eq. (46). The product of two four dimensional representations 
(ai, a2, a 3 , a 4 ) ~ 4_i and (61, b 2 , b 3 , b^) ~ 4^ contains the singlet 



(ab) : 



(ai&i + a 2 6 2 + «3^3 + 0464) 



(47) 



and the triplets: 



\ 



-a 4 6i + a 3 b 2 - a 2 b 3 + ai& 4 
-a 3 h - a 4 b 2 + aib 3 + a 2 b A J (06)3, 
a 2 h - a\b 2 - a 4 & 3 + a 3 6 4 

a\b\ - a 2 b 2 - a 3 b 3 + a 4 6 4 
-aih + a 2 b 2 - a 3 b 3 + a 4 b A ] (a&) 3/ 
-a\b\ - a 2 b 2 + a 3 b 3 + a 4 6 4 



' — a 4 6i — a 3 b 2 + a 2 6 3 + ai& 4 
a 3 6i — a 4 &2 — aib 3 + a 2 bi 
1 -a 2 6i + ai6 2 - a 4 6 3 + a 3 6 4 



^ a 4 6i + a 3 6 2 + a 2 b 3 + ai& 4 

a 3 6i + a 4 6 2 + ai^3 + «2&4 
y a 2 6i + ai6 2 + a 4 6 3 + a 3 6 4 

^ a 4 6i - a 3 5 2 - a 2 6 3 + ai& 4 ^ 

-a 3 6i + a 4 6 2 - ai^3 + «2&4 
y -a 2 b\ - a\b 2 + a 4 6 3 + a 3 6 4 y 

(48) 
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